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Abstract 

We extend a recent chiral approach to nuclear matter of Lutz et al. [Phys. Lett. 
B474 (2000) 7] by calculating the underlying (complex-valued) single-particle poten- 
tial U{p,kf) + iW{p,kf). The potential for a nucleon at the bottom of the Fermi-sea, 
U{0,kfQ) = — 20.0MeV, comes out as much too weakly attractive in this approach. Even 
more seriously, the total single-particle energy does not rise monotonically with the nu- 
cleon momentum p, implying a negative effective nucleon mass at the Fermi-surface. Also, 
the imaginary single-particle potential, W{0, kfo) = 51.1 MeV, is too large. More realistic 
single-particle properties together with a good nuclear matter equation of state can be 
obtained if the short range contributions of non-pionic origin are treated in mean-field ap- 
proximation (i.e. if they are not further iterated with Ivr-exchange) . We also consider the 
equation of state of pure neutron matter En{kn) and the asymmetry energy A{kf) in that 
approach. The downward bending of these quantities above nuclear matter saturation 
density seems to be a generic feature of perturbative chiral pion-nucleon dynamics. 

PACS: 12.38.Bx, 21.65.+f 
Keywords: Nuclear matter equation of state; Complex single-particle potential in symmetric 
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1 Introduction 

The present status of the nuclear matter problem is that a quantitatively successful description 
can be achieved, using advanced many-body techniques ||l|, in a non-relativistic framework 
when invoking an adjustable three-body force. Alternative relativistic mean-field approaches, 
including non-linear terms with adjustable parameters or explicitly density- dependent point 
couplings, are also widely used for the calculation of nuclear matter properties and finite nuclei 

i- 

In recent years a novel approach to the nuclear matter problem based on effective field 
theory (in particular chiral perturbation theory) has emerged [0, ^ ^. The key element there 
is a separation of long- and short- distance dynamics and an ordering scheme in powers of small 
momenta. At nuclear matter saturation density the Fermi-momentum kjQ and the pion mass 
are comparable scales {kfo ^ 2m^), and therefore pions must be included as explicit degrees of 
freedom in the description of the nuclear many-body dynamics. The contributions to the energy 

1 Work supported in part by BMBF, GSI and DFG. 



per particle of isospin-symmetric nuclear matter E{kf) originating from chiral pion-nucleon 
dynamics have been calculated up to three-loop order in refs. 0, ^ . Both calculations include 
the Ivr-exchange Fock-diagram and the iterated Ivr-exchange Hartree- and Fock-diagrams. In 
ref.0 irreducible 27r-exchange is also taken into account and a momentum cut-off A is used 
to regularize the few divergent parts associated with chiral 27r-exchange. The resulting cut-off 
dependent contribution to E{kf) is completely equivalent to that of a zero-range NN-contact 
interaction (see eq.(15) in ref.0]). At that point the work of Lutz et al. |^ deviates and it 
follows a different strategy. Two zero-range NN-contact interactions (acting in ^5*1 and ^5*0 
NN-states) proportional to the parameters 5^0 + 5^1/4 and gi+g\/4: are introduced (see eq.(4) in 
ref.p]). The components proportional to g\/4: cancel the zero-range contribution generated by 
the Ivr-exchange Fock-diagram. The other components proportional to Qq and gi are understood 
to subsume all non-perturbative short-range NN-dynamics relevant at densities around nuclear 
matter saturation density po- In order to be consistent with this interpretation the NN-contact 
vertices proportional to (70,1 ^-re allowed to occur only in first order. Furthermore, according to 
ref.0 pions can be treated perturbatively (at least) in the ^5*0 partial-wave of NN-scattering 
if the zero-range pieces they generate are removed order by order. Therefore, the NN-contact 
vertex proportional to g\/4: occurs also in higher orders (see Fig.l in ref. p| which includes 
diagrams with "filled circle" and "open circle" vertices). 

Despite their differences in the treatment of the effective short-range NN-dynamics both 
approaches p, |[ are able to reproduce correctly the empirical nuclear matter properties (sat- 
uration density pq, binding energy per particle —Eikf^) and compressibility K) by adjusting 
only one parameter, either the coupling (/q + (71 ~ 3.23 or the cut-off A ~ 0.65 GeV. Note that 
in dimensional regularization all diagrams evaluated in ref.[|l are finite. In the chiral approach 
of the Munich group 0, Q the asymmetry energy A{kf), the energy per particle of pure neutron 
matter En{kn) as well as the (complex) single-particle potential f/(p, kf) +i W{p, kj) below the 
Fermi-surface (p < kf) have been calculated. Good results (in particular for the asymmetry en- 
ergy, A{kfo) = 33.8 MeV, and the depth of the single-particle potential, U (0, fc/o) = —53.2 MeV) 
have been obtained with the single cut-off scale A ~ 0.65 GeV adjusted to the binding energy 
per particle —Eikfo) = 15.3 MeV. Moreover, when extended to finite temperatures this ap- 
proach reproduces the liquid-gas phase transition of isospin-symmetric nuclear, however, with 
a too high value of the critical temperature Tc = 25.5 MeV. 

It is the purpose of this work to investigate in the approach of Lutz et al. ^ the single-particle 
potential U{p, kf)+i W{p, kf) in isospin-symmetric nuclear matter as well as the neutron matter 
equation of state En{kn) and the asymmetry energy A{kf). One of our major conclusions will 
be that any strong short-range NN-dynamics of non-pionic origin should be kept at the mean- 
field level. It should not be further iterated with pion-exchange in contrast to the prescription 
of power-counting rules. 



2 Nuclear matter equation of state 

Let us first reconsider the equation of state of isospin-symmetric nuclear matter as it follows 
from the calculation of ref.0. Even though all contributions to the energy per particle E{kf) 
have been given explicitly in ref . ||^ we prefer to write down again the extra terms generated by 
the NN-contact interactions proportional to g^^i + g\/4 (using a more compact notation). The 
first diagram in Fig. 1 gives rise to a contribution to energy per particle of the form: 

s,, , (7 + 1)9^*:? „, 
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Figure 1: Additional in-medium diagrams generated by the NN-contact interactions introduced 
in ref.P]. The two NN-contact interactions proportional to 7 -|- 1 and 7™ -|- 1 are symbolized by 
the filled square vertex. The last diagram is to be understood such that quadratic terms (such 
as 7^,77n and 7^) are omitted. 



where we have introduced (for notational convenience) the coefficient 7 by the relation (7 + 
l)g\/2 = go + gi + g\/2. In the second and third diagram in Fig. 1 the contact-interaction 
proportional to 7 + 1 is iterated with Ivr-exchange or with itself (dropping the 7^-contribution). 
Putting a medium-insertionQ at each of two nucleon propagators with equal orientation one 
gets: 
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with the abbreviation u = kf/niT^. One observes that eq.(2) receives no contribution from the 
third diagram in Fig. 1 since dl 1 is set to zero in dimensional regularization. The second 
and third diagram in Fig. 1 with three medium-insertions give rise to the following contribution 
to the energy per particle: 
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with the auxiliary functions H = \n{u + xy) — \n{u — xy) and s = xy + a/m^ — + x^y^. In 
the chiral limit = only the contribution coming from the last term, — in the second 
square bracket survives. The corresponding double integral dx x'^ J^^ dy . . . has the value 
2-u^(ln4 — 11)/105. The expansion of the energy per particle up to order 0{kj) is completed 
by adding to the terms eqs. (1,2,3) the contributions from the (relativistically improved) kinetic 
energy, from Ivr-exchange and from iterated Ivr-exchange written down in eqs. (5-11) of ref. [Q. In 
case of the Ivr-exchange contribution (eq.(6) in ref. [^]) we neglect of course the small relativistic 
1/M^-correction of order 0{kj). 

Now, we have to fix parameters. The pion decay constant /^r = 92.4 MeV and the nucleon 
mass M = 939 MeV are well-known. As in ref.0] we choose the value gA = 1-3. This corre- 
sponds via the Goldberger-Treiman relation to a vriViV-coupling constant of g-^N = gAM/ = 
13.2 which agrees with present empirical determinations of g-j^N from TrA^-dispersion relation 
analyses [§. We set m,r = 135 MeV (the neutral pion mass) since this is closest to the expected 
value of the pion mass in the absence of isospin-breaking and electromagnetic effects. 

The dashed line in Fig. 2 shows the equation of state of isospin-symmetric nuclear matter in 
the approach of ref. 0] using the abovementioned input parameters. The coefficient 7 = 4.086 
has been adjusted such that the minimum of the saturation curve lies at E{kfQ) = —15.3 MeV 
0. The predicted equilibrium density po = 0.138 fm~^ (corresponding to a Fermi-momentum of 

^This is a technical notation for the difference between the in-medium and vacuum nucleon propagator. For 
further details, see section 2 in ref. B. 
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Figure 2: The energy per particle of isospin-symmetric nuclear matter E{kf) versus the nucleon 
density p = 2kj/3n'^. The dashed line corresponds to the approach of ref. P]. The full line results 
if the NN-contact interaction is treated in mean-field approximation. In each case the coefficient 
7 is adjusted such that the saturation minimum lies at Ei^kjo) = — 15.3MeV. 



kfo = 250.1 MeV) is somewhat too low. The same holds for the nuclear matter compressibility 
K = kjoE"{kfo) = 202 MeV. Of course, if we use the input parameters of ref.|l (/^ = 93 MeV, 
gA = 1.26, m,r = 140 MeV and go + di = 3.23 corresponding to 7 = 4.07) we exactly reproduce 
the numerical results of that work. We emphasize that the different treatment of the two 
components of the NN-contact interaction is essential in order to get (realistic) nuclear binding 
and saturation in the framework of ref . . If both components were treated on equal footing 
in first order (technically this is realized by deleting the contribution coming from the third 
diagram in Fig. 1) the energy per particle E{kf) would not even develop a minimum. 



3 Real single-particle potential 

Next, we turn to the real part of the single-particle potential U{p, kf) below the Fermi-surface 
{p < kf) in the framework of ref.p|. As outlined in ref. the contributions to U{p, kf) can be 
classified as two-body and three-body potentials. From the first diagram in Fig. 1 one gets a 
contribution to the two-body potential of the form: 
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Figure 3: The lower curve shows the real part of the single-particle potential U{p,kfQ) at 
saturation density /c/o = 250.1 MeV in the approach of Lutz et al. f^. The upper curve includes 
in addition the relativistically improved kinetic energy Tkin(p) = p^/2M — p^/8M^. 



which is just twice its contribution to the energy per particle (see eq.(l)). From the second 
diagram in Fig. 1 one derives a contribution to the two-body potential of the form: 
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with the abbreviation x = p/rriT^. The second and third diagram in Fig. 1 give each rise to three 
different contributions to the three-body potential. Altogether, they read: 
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with the auxiliary function a = iy+^Ju^ ~ + ^'^v'^- The real single-particle potential U {p, kf) 
is completed by adding to the terms eqs. (4,5,6) the contributions from Ivr-exchange and iterated 
Ivr-exchange written down in eqs. (8-13) of ref.|p. Again, the (higher order) relativistic 
correction to Ivr-exchange (see eq.(8) in ref.]^) is neglected for reasons of consistency. 
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Figure 4: The lower curve shows the real part of the single-particle potential U{p,kfQ) at 
saturation density /c/o = 270.3 MeV in a mean- field treatment of the NN-contact interaction. 



The lower curve in Fig. 3 shows the momentum dependence of the real single-particle 
potential U{p,kfQ) at saturation density kfo = 250.1 MeV as it arises in the framework of 
Lutz et al. 0. The predicted potential depth for a nucleon at the bottom of the Fermi-sea 
is only f/(0, fc/o) = —20.0 MeV. In magnitude this is much smaller than the typical depth 
Uq — —53 MeV of the empirical optical model potential [0 or the nuclear shell model po- 
tential |]ll| . The upper curve in Fig. 3 includes the (relativistically improved) single-nucleon 
kinetic energy Tkin(p) = /2M — p'^/8M^. As required by the Hugenholtz- van- Hove theorem 
[1T2| this curve ends at the Fermi-surface p = kfo with the value E{kfQ) = —15.3 MeV. A further 
important check is provided by the sum rule for the two- and three-body potentials f/2,3(p, kf) 
written down in eq.(5) of ref.|^. It holds with very high numerical accuracy in the present 
calculation. 

The momentum dependence of the two (dashed) curves in Fig. 3 is completely unrealistic. 
Most seriously, the total single-particle energy Tkin(p) + U{p,kfo) (upper curve) does not rise 
monotonically with the nucleon momentum p, but instead it starts to bend downward above 
p ^ 190 MeV. This implies a negative effective nucleon mass at the Fermi-surface, M*{kfo) ~ 
—3.5 M, and a negative density of states with dramatic consequences for the finite temperature 
behavior of nuclear matter. Because of such pathological features of the underlying single- 
particle potential the scheme of Lutz et al. [Q has to be rejected in its present form. 

The overly strong momentum dependence of U{p,kfQ) comes from the second and third 
diagram in Fig. 1 in which the NN-contact interaction proportional to the large coefficient 7-I- 1 
is further iterated. We propose to drop these three-loop diagrams and to keep the NN-contact 
interaction (of unspecified dynamical origin) at the mean-field level. The resulting equation of 
state obtained by leaving out the contributions eqs.(2,3) and adjusting 7 = 6.198 is shown by 
the full line in Fig. 2. The predicted saturation density is now po = 0.174 fm~'^ (corresponding 
to a Fermi-momentum of kfo = 270.3 MeV) and the nuclear matter compressibility has the 



value K = 253 MeV. Note that the scheme of ref . modified by a mean-field treatment of the 
NN-contact interaction becomes equivalent to the truncation at fourth order in small momenta 
of our previous work |^ ^ after the identification of parameters, 7 + 1 = 10g\AM / {4it f-j^Y , 
with A denoting the cut-off scale. 

The lower full curve in Fig. 4 shows the momentum dependence of the real single-particle 
potential at saturation density kfo = 270.3 MeV which results in a mean-field approximation of 
the NN-contact interaction (by leaving out the contributions eqs.(5,6)). The predicted potential 
depth U{0,kfo) = —54.8 MeV is in good agreement with that of optical model |10] or nuclear 
shell model potentials |Tl|]. Most importantly, the total single-particle energy Tkin(p) + f/(p, fc/o) 



(upper curve) grows now monotonically with the nucleon momentum p, as it should. The up- 
and downward bending of the lower full curve in Fig. 4 is however still too strong. The negative 
slope of f/(p, fc/o) at the Fermi-surface p = fc/o leads to a too large effective nucleon mass 
M*{kfQ) ^ 2.9 M which reflects itself in a too high critical temperature ~ 25 MeV of the 
liquid-gas phase transition . More elaborate calculation of nuclear matter in effective (chiral) 
fleld theory are necessary in order to cure this problem of the too large effective nucleon mass 

M*(fc;o). 



4 Imaginary single-particle potential 



In this section, we discuss the imaginary part of the single-particle potential W{p, kf) for p < kf 
as it arises in the scheme of Lutz et al. 0. This quantity determines the half-width of nucleon- 
hole states in the Fermi-sea. As outlined in ref.[^ the contributions to W{p, kf) can be classifled 
as two-body, three-body and four-body terms. From the second and third diagrams in Fig. 1 
one derives a two-body term of the form: 
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The associated three-body term reads: 
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and the four-body term is given by the expression: 
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Figure 5: The imaginary part of the single-particle potential W{p,kfQ) at saturation density 
versus the nucleon momentum p. The dashed line corresponds to the approach of Lutz et 
al. 1^ and the full line shows the result obtained in a mean-field treatment of the NN-contact 
interaction. 



The additional contributions from the iterated Ivr-exchange Hartree- and Fock-diagram are 
collected in eqs. (20-25) of ref.p. The total imaginary single-particle potential evaluated at 
zero nucleon momentum {p = 0) can even be written as a closed form expression: 



W{0,kj) 



u 



+ 
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where Li2(— ci~^) = + a)~^\ii( denotes the conventional dilogarithmic function. 

The dashed line in Fig. 5 shows the momentum dependence of the imaginary single-particle 
potential W{p, kfo) at saturation density fc/o = 250.1 MeV as it arises in the approach of ref. p. 
The predicted value W{0,kfQ) = 51.1 MeV lies outside the range 20 — 40 MeV obtained in 
calculations based on (semi)-realistic NN- forces [jl3|, . The full line in Fig. 5 corresponds 
to a mean-field approximation of the NN-contact interaction. Up to a slight change in the 
equilibrium Fermi-momentum kfo = 270.3 MeV the full curve in Fig. 5 agrees with the one 
shown in Fig. 4 of ref. 0. The considerably reduced value W{0, fc/o) = 28.4 MeV indicates the 
large contribution of the iterated diagrams in Fig. 1 to the imaginary single-particle potential 
W{p, kf). Note that both curves in Fig. 5 vanish quadratically near the Fermi-surface as required 
by Luttinger's theorem [|15]. As further check on our calculation we verified the zero sum 
rule: J^^^ dpp'^[6W2{p,kf) + 4PV3(p, fcj) + SW^^p, kf)] = 0, for the two-, three- and four-body 
components ^2,3,4(^5^/) written in eqs. (7,8,9). 
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Figure 6: The energy per particle of pure neutron matter En{kn) versus the neutron density 
Pn = ^n/3vr^. The dashed hne corresponds to the approach of Lutz et al. and the full 
line shows the result obtained in mean-field approximation of the nn-contact interaction. The 
dotted line stems from the many-body calculation of the Urbana group ||16 



5 Neutron matter 

In this section we discuss the equation of state of pure neutron matter. In the scheme of Lutz 
et al. [§] the energy per particle of pure neutron matter En{kn) depends exclusively on the 
coefficient gi parameterizing the short-range NN-interaction in the channel with total isospin 
1 = 1. There is no need to write down explicitly the contributions of the diagrams in Fig. 1 to 
Enikn). These expressions are easily obtained from eqs. (1,2,3) by replacing kf by the neutron 
Fermi-momentum by replacing the coefficient 7 by a new one 7^, and by multiplying the 
formulas with a relative isospin factor 1/3. The relation (7^ -|- l)g\/4: = gi + fi'1/4 defines 
this new coefficient 7„. The additional contributions to En{kn) from the kinetic energy, Ivr- 
exchange and iterated Ivr-exchange are written down in eqs. (32-37) of ref.[^ (neglecting again 
the relativistic 1/M^-correction to Ivr-exchange). 

The dashed line in Fig. 6 shows the energy per particle of pure neutron matter En{kn) versus 
the neutron density pn = kfj^/Sn^ as it arises in the approach of ref.0. The coefficient 7„ = 0.055 
has been adjusted to the empirical value of the asymmetry energy A{kfQ = 250.1 MeV) = 
33.2 MeV (see next section). The downward bending of the dashed curve in Fig. 6 above p„ > 
0.15 fm~^ is even stronger than in our previous work (see Fig. 8 therein). This property can 
be understood by taking the chiral limit (m^r — > 0) of the calculated neutron matter equation 
of state and considering the coefficient /3„ in front of the term k^/M^. In the approach of Lutz 
et al. one has: 

= + 17+ 16 In 2) - 1 = -1.23, (11) 

which is 2.2 times the negative value of /?„ found in ref. |Q. The full line in Fig. 6 shows the 



equation of state of pure neutron matter obtained in mean-field approximation of the nn- 
contact interaction proportional to 7„ + 1 after adjusting 7„ = 0.788 to the empirical value 
of the asymmetry energy A{kfQ = 270.3 MeV) = 33.2 MeV. The downward bending of the 
full curve in Fig. 6 is weaker and it sets in at somewhat higher densities p„ > 0.2 fm~'^. The 
dotted line in Fig. 6 stems from the many-body calculation of the Urbana group |jl6|. This 



curve should be considered as a representative of the host of existing realistic neutron matter 
calculations which scatter around it. The systematic deviations observed in Fig. 6 indicate that 
the neutron matter equation of state of ref. |^ cannot be improved by a different treatment of 
the short-range NN-dynamics alone. The downward bending of En{kn) above p„ > 0.2 fm~^ 
seems to be a generic feature of perturbative chiral vrA^-dynamics truncated at three-loop order. 



6 Asymmetry energy 

Finally, we turn to the density dependent asymmetry energy A{kf) in the approach of ref. [Q]. 
The asymmetry energy is generally defined by the expansion of the energy per particle of isospin- 
asymmetric nuclear matter (described by different proton and neutron Fermi momenta kp^n = 
kf{l =F Sy^^) around the symmetry line: E^{kp, = E{kf) + 5^ A{kf) + 0{5'^). Evaluation of 
the first diagram in Fig. 1 leads to the following contribution to the asymmetry energy: 
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with the coefficients '~i = '^{gQ + gi) / g\ and 7„ = Agi/g\ in the notation of ref. 0. Putting a 
medium-insertion at each of two nucleon propagators with equal orientation one gets from the 
second and third diagram in Fig. 1: 



A{kf) 



3(87r)3/4 



^ , 2(7 + l)u + 8(27 -7„ + l)u^ arctan 2m 



+ 



(27„ - 67 - A)u - 



7 + 1 
2u 



ln(l + W 



(13) 



The same diagrams with three medium-insertions give rise to the following contribution to the 
asymmetry energy: 
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with s' = uds/du and s" = u'd'^s/du' denoting partial derivatives. In the chiral limit 
m^r = only the terms in the first line of eq.(14) survive. The corresponding double inte- 
gral dx x'^ J^^ dy . . . has the value 4M^(ln4 — 1)/15. The asymmetry energy is completed 
by adding to the terms eqs.(12,13,14) the contributions from the kinetic energy, (static) Itt- 
exchange and iterated Ivr-exchange written down in eqs. (20-26) of ref. [§]. 




Figure 7: The asymmetry energy A{kf) versus the nucleon density p = 2/cj/37r^. The dashed 
hne corresponds to the approach of Lutz et ah and the full line shows the result obtained 
in mean-field approximation of the NN-contact interactions. The parameter 7„ is in each case 
adjusted to the (empirical) value A(/c/o) = 33.2 MeV 0. 

The dashed line in Fig. 7 shows the density dependence of the asymmetry energy A{kf) in 
the approach of Lutz et al. [Q] with the coefficient 7„ = 0.055 adjusted (at fixed 7 = 4.086) to 
the empirical value A[kfo = 250.1 MeV) = 33.2 MeV 0. The full line in Fig. 7 corresponds to 
the result obtained in mean-field approximation of the NN-contact interaction by dropping the 
contributions eqs. (13,14). In that case the empirical value A{kfQ = 270.3 MeV) = 33.2 MeV [H 
is reproduced by tuning (at fixed 7 = 6.198) the coefficient 7„ to the value 7„ = 0.788. Both 
curves in Fig. 7 behave rather similarly. In each case the asymmetry A{kf ) reaches it maximum 
close to the respective saturation density po and then it starts to bend downward. Since the 
same (unusual) feature has also been observed in ref . it seems to be generic for perturbative 
chiral vriV-dynamics truncated at three-loop order. 

7 Concluding remarks 

In this work have we continued and extended the chiral approach to nuclear matter of Lutz 
et al. by calculating the underlying single-particle potential. The potential for a nucleon 
at the bottom of the Fermi-sea U{0,kfo) = —20.0 MeV is not deep enough. Most seriously, 
the total single-particle energy Tkin(p) + U{p,kfQ) does not grow monotonically with the nu- 
cleon momentum p. The thereof implied negative effective nucleon mass at the Fermi-surface 
M*{kfQ) — 3.5M and the negative density of states will ruin the behavior of nuclear mat- 
ter at finite temperatures. The half-width of nucleon-holes at the bottom of the Fermi-sea 
W{0,kfo) = 51.1 MeV comes also out too large in that approach. A good nuclear matter equa- 
tion of state and better (but still not yet optimal) single-particle properties can be obtained if 



the NN-contact interaction (proportional to the coefficient go + Qi +fl'l/2) is kept at the mean- 
field level and not further iterated. The energy per particle of pure neutron matter and 
the asymmetry energy A{kf) depend on a second parameter gi in the scheme of ref. [0. Their 
density dependence is similar to the results of the one-parameter calculation in ref. 0. The 
downward bending of En{kn) and A{kf) above saturation density po (less pronounced if the 
NN-contact interaction is kept at mean-field level) seems to be generic for perturbative chiral 
TrA^-dynamics. More elaborate calculations of nuclear matter in effective (chiral) field theory 
which fulfill all (semi)-empirical constraints are necessary. 
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